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Abstract 

Firstly, a sequential symmetry adaptation procedure is derived for semidirect product groups. 



,-S^ . Then, this sequential symmetry adaptation procedure is used in the development of new method 

named Monomer Basis Representation (MBR) for calculating the vibration-rotation-tunneling 

(VRT) spectra of molecular clusters. The method is based on generation of optimized bases for 
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each monomer in the cluster as a linear combination of some primitive basis functions and then 
using the sequential symmetry adaptation procedure for generating a small symmetry adapted ba- 
sis for the solution of the full problem. It is seen that given an optimized basis for each monomer 
the application of the sequential symmetry adaptation procedure leads to a generalized eigenvalue 

problem instead of a standard eigenvalue problem if the procedure is used as it is. In this paper, 
> ■ 
^C MBR method will be developed as a solution of that problem such that it leads to generation of an 



orthogonal optimized basis for the cluster being studied regardless of the nature of the primitive 



bases that are used in the generation of optimized bases of the monomers. 
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I. INTRODUCTION 

Obtaining the bulk properties of condensed phases from microscopic considerations is 
a long standing dream. For that purpose, it is essential to develop an accurate potential 
surface. 

The total potential energy function V(f[, r~2, . . . , rjv) of a system of N identical particles 
can be expanded as sum of n-body potentials: 

N N 

V(fl, f 2l ■ ■ ■ , fff) = Y^ ^( r *i> 'w) + Yl ^3(^1,^2,^3) 

il<l2 il<«2<«3 

N 
"T • • • ~r / Vnviu ^2, ^3, • • ■ , '"jjv/ - V / 

U<i2<«3,---i*iV 

It has been assumed in the past that this series converges rapidly and in most systems the 
pairwise additive approximation, i.e. 

N 

V(f{, f 2 ,...,fN)= Y Viirti, rQ, (2) 

h<i2 

ias been shown to be qualitatively valid. However, this is not the case for many systems 
lH4|. For example, the higher order terms in equation (TjQ) may add up to 25% of the 
interaction energy of bulk water [5J, most of which results from three-body effects. It has 
been found that 6| the pair interaction energy represents 83 — 86% of the total interaction 
energy of (H 2 0)3 ~ 75% of (H 2 0)± and only w 68% of {H 2 0)^. Therefore, it is evident that 
many-body terms are far from negligible. 

Before the studies of van der Waals molecules, most experiments sensitive to many-body 
effects were performed on macroscopic systems with the result that only the total many 
body effects were obtained. Then, condensed phase properties have been used to determine 
"effective" pair potentials |7|] which include all of the terms in the many-body expansion in 
an average manner, at least over some range of temperature. 

On the other hand, studies of small clusters proved to be ideal for studying many-body 
terms, because each term in the many-body expansion can be studied uniquely. For example, 
dimers are ideal for studying two-body potentials, since there are no many-body terms in 
their potential surfaces. Similarly, once the potential surface of the dimer of a substance 
is known accurately, the trimer of that substance becomes the model system for studying 



three-body terms in its potential surface, because three-body terms are the only many- 
body terms in the potential surface of a trimer. Thus, in principle it is possible to develop 
accurate potential surfaces by starting from dimers and deriving the many-body terms by 
studying bigger and bigger clusters. However, it is very hard to accomplish in practice. In 
order to develop a potential surface or in order to test the accuracy of an available potential 
surface, it is necessary to make quantum mechanical calculations. In quantum mechanical 
spectra calculations, the eigenstates of the Hamiltonian operator are usually obtained by 
diagonalizing the matrix representing the Hamiltonian operator in a finite basis and it is 
well known that the size of the basis required for the calculations grows exponentially with 
the number of degrees of freedom. As a result of that, the quantum mechanical calculations 
becomes harder as the size of the problem gets bigger. 

Considering inter-molecular vibrations of molecular clusters (consisting of nonlinear 
monomers), dimers are six dimensional systems and trimers are twelve dimensional systems. 
If it is possible to obtain accurate results with N basis functions per degree of freedom, then 
a dimer calculation requires iV 6 basis functions and a trimer calculation requires N 12 basis 
functions. Even if it would be possible to obtain good results with a small number of basis 
functions per degree of freedom such as N = 10, the size of a basis required for a trimer 
calculation would be a million times bigger than the size of a basis required for a dimer 
calculation. It is this exponential scaling of the number of the basis functions which makes 
the calculations harder and harder as the size of the problem gets bigger, and obviously it 
is necessary to develop methods for reducing the size of the bases. 

One way of reducing the sizes of bases is to make symmetry adaptation of basis functions 
and solve for each symmetry separately. However, the well-known method of symmetry 
adaptation is not very helpful. Although, the sizes of bases grow exponentially, orders of 
molecular symmetry groups [8j grow linearly. For example, order of the molecular symmetry 
group of water dimer is 16 [9J and the order of the molecular symmetry group of water trimer 
is 48 [id]. 

A much more efficient way of reducing the computational cost is to use bases that are 
optimized for the particular problem at hand instead of primitive bases. The inefficiency of 
primitive bases results from the fact that they do not know anything about the potential 
surface of the system. The optimized bases that know about the underlying potential surface 
can be generated as linear combinations of some primitive basis functions by taking a model 



potential surface that resembles to the actual potential surface as much as possible and 
solving for its eigenstates. This obviously makes it necessary to divide the problem into 
smaller parts since trying to find an optimized basis for the full problem is as difficult as 
solving it. In the case of molecular clusters an obvious way of dividing the problem into 
smaller parts is to consider each monomer separately. In the rest of the paper, the main 
discussion will be about how optimized bases for each monomer can be generated and how 
they can be combined for the solution of the full problem. 

In section [Til a sequential symmetry adaptation procedure will be derived. By finding the 
relations between the projection operators of the irreducible representations of the molecular 
symmetry group of the cluster and the projection operators of the irreducible representations 
of its subgroups, symmetries of monomer basis functions will be related to the symmetries 
of eigenstates of the cluster. This symmetry adaptation procedure will be used in section ITTT1 
in which generation of optimized monomer bases is discussed. It will be seen that generation 
of optimized bases creates its own problems related with symmetry adaptation. In order 
to guarantee generation of an efficient orthonormal basis, it will be necessary to modify 
the sequential symmetry adaptation procedure developed in section [Til In order to illustrate 
how the method can be used for generating small symmetry adapted bases, two qualitative 
models will be given in section [V] one for water dimer and one for water trimer. The method 
will be applied to water dimer in a forthcoming paper in order to illustrate its application 
and to show its validity. 

II. SEQUENTIAL SYMMETRY ADAPTATION 

An analysis of the structure of the molecular symmetry groups shows that the molecular 
symmetry group of a molecular cluster consisting of n non-reacting monomers can be written 



in terms of its subgroups as H] 



G(MS) = ((G kl ®G k2 ®...® G kn ) ©d) ® e. (3) 

In the equation above, the groups G ki with i = 1, . . . ,n are the pure permutation subgroups 
of the monomers that have orders kf, the group Gi which is of the order of I is the subgroup 
containing the operations permuting the identical monomers; and the group e is the inversion 
subgroup that contains the identity element E and the inversion element E*. 



In equation ([3]), ® denotes a direct product multiplication and (g) denotes a semi-direct 
product multiplication. The difference between a direct product and a semi-direct product 
multiplication is that in a direct product multiplication both of the subgroups are invariant 
subgroups of the product group while in a semi-direct product multiplication only one of 
them is an invariant subgroup of the product group. Since the operations permuting the 
identical monomers bring in non-commutation, presence of a semi-direct product multipli- 
cation is inevitable. In this paper, it will be assumed that the group Gi is cyclic and has the 
order / = n m . In the case of dimers the order of the group which contains the permutations 
of identical monomers has to be 2 anyway (since n m \ = n m for n m = 2), so that there is no 
assumption. However, in the case of bigger clusters, physical meaning of this assumption is 
that the cluster is assumed to have a rigid cyclic skeleton. The reason for this assumption 
will be clear in section IIIIB1 It should be noted that equation fl3]) assumes that the inversion 
operation is a feasible symmetry operation of the molecular symmetry group which is not 
true in general. However, this corresponds to a more general case. The method developed 
in this paper is applicable to any cluster whose molecular symmetry group does not include 
inversion operation, too. 

Symmetry adaptation of basis functions to an irreducible representation T of a group 
G can be done by application of the projection operator of that irreducible representation 



which is given by 12|, |l3] 



^ r = ^£Am, (4) 

|G| a eG 
where g are the elements of the group G, dr is the dimension of the irreducible representation 
T, \G\ is the order of the group G, O g is the operator representing g, and x r [d] is the character 
of g in the irreducible representation T. As shown in the appendix, for a group G that can 
be written as a semi-direct product of two of its subgroups H and K, and satisfying the 
condition given in equation flHj), the projection operator of an irreducible representation 
can be decomposed into product of two terms such that 



r G 



wX^yoMwX^A' (5) 

s ' ' heH / \ ' ' k£K J 



where the definitions of the terms are similar to those of equation (j4]). Please note that the 
condition given in equation ( pH| is always satisfied for any direct product group. Therefore, 
equation ([5]) applies to all of the group multiplications in equation fl3]). Consequently, it 



follows that the symmetry adaptation of basis functions can be done in n + 2 steps sequen- 
tially Furthermore, the characters in equation (JSJ) can be decomposed to their irreducible 
components in groups H and K. Thus, if the character in the first parentheses in equation 
(E]), can be expressed as a\Y\ © a 2 T 2 © . . . where a, are nonnegative integers and Tj are the 
irreducible representations of the group H, then from the definition of projection operators 
it follows that this parentheses can be expressed in terms of the projection operators of the 
irreducible representations of the group H as 

ij2x T [h]*d h = a -±pZ + fP2 + ..., (6) 



\H\ ^^ L J d x n d 

1 ' heH 



il\ 



where P H % is the projection operator and di are the dimension of the irreducible represen- 
tation Tj. A similar expression can also be written for the second parentheses in equation 
(JSJ) in terms of the projection operators of the irreducible representations of the group K. 
Obviously only the terms for which dj is nonzero will contribute to the sum. Consequently, 
equation (JBJ), or equally finding the coefficients a\ is useful for determining which symmetries 
of the monomer basis functions can be symmetry adapted to an irreducible representation 
of the product group. 

III. MONOMER BASIS REPRESENTATION METHOD 

If the symmetry adaptation procedure given in previous section is used with primitive 
bases, it can not provide any optimization more than what one can achieve with the di- 
rect symmetry adaptation by using equation (J4]). On the other hand, sequential symmetry 
adaptation procedure combined with the physically meaningful partitioning of the molecular 
symmetry groups given in equation ((3]) makes it possible to devise algorithms for obtaining 
symmetry adapted optimized bases. 

On the other hand, generation of optimized bases creates its own problems related with 
symmetry adaptation. Primitive basis functions (plane waves, spherical harmonics, Wigner 
rotation functions . . . ) are the solutions of Hamiltonians corresponding to motions of free 
particles or free bodies. Since the kinetic energy is always absolutely symmetric, a free 
particle Hamiltonian has absolute symmetry, too. Consequently, the basis functions that 
are obtained as solutions of that Hamiltonian also have absolute symmetry. As a result 
of that application of a symmetry operation to a primitive basis function always results in 

6 



another function in the same basis other than a possible phase factor. Thus the symmetry 
adapted basis functions can be obtained as linear combinations of primitive basis functions. 

The case of optimized basis functions is different. An optimized basis function should 
know about the particular problem at hand so that the Hamiltonian, that the optimized basis 
functions are solutions of, should include a potential energy function. Since the potential 
surfaces do not have absolute symmetry, a Hamiltonian including a potential energy function 
can not be absolutely symmetric, either. This will restrict the symmetries of optimized basis 
functions that are obtained as solutions of the Hamiltonian. This means that application of 
a symmetry operation to an optimized basis function will not necessarily result in another 
basis function in the same optimized basis. As a result of that solving a problem might 
be impossible when optimized basis functions are used unless a special care is taken to 
ensure that the physically meaningful solutions of the Hamiltonian can be obtained as linear 
combinations of the optimized basis functions. 

A discussion of how such optimized monomer bases can be found and how they can be 
combined for the solution of the full problem will be made in the following subsections. 
The procedure for obtaining the results can be divided into four steps. These steps are: 
generation of a basis for one of the monomers, generation of bases for other monomers, 
combining monomer bases to generate a basis for the full problem, and the solution of the 
full problem. 

Before starting to develop the method, it should be noted that a basis function related 
with a monomer is a function which describes the orientation of the monomer in the cluster 
and a function related with inter-monomer coordinates is a function which describes the 
orientation of the monomers with respect to each other. The monomers are assumed to be 
rigid bodies so that intra-monomer degrees of freedom are not considered. In the discussion, 
language of Permutation Inversion (PI) group theory [8|, [l2|, ll4j will be used since it provides 
the most natural way of handling the symmetries in molecular systems. 

A. Generation of A Monomer Basis 

In order to have a basis which is fully symmetry adapted, or can be fully symmetry 
adapted, to the molecular symmetry group of the cluster, there should be a basis for each 
monomer, which is symmetry adapted to the irreducible representations of the pure per- 



mutation group of the monomer Gk- This can be done by just taking a primitive basis 
for each monomer, and then symmetry adapting it to the irreducible representations of the 
permutation group, Gk- However, in order to have an efficient basis for the full problem, 
the size of the monomer basis should be small; and just taking a small number of primitive 
basis functions for each monomer will not be a good choice obviously. Therefore, instead of 
using a primitive basis for each monomer, which does not know anything about the problem 
at hand, it is better to have a basis which is optimized for the particular problem being 
studied. 

An optimized basis for a monomer can be generated by taking a model Hamiltonian for 
that monomer and then solving for the eigenstates of the model Hamiltonian with a basis 
which has the required symmetry properties. Then, a small number of the eigenstates of 
the model Hamiltonian can be taken as an optimized basis for that monomer. The model 
Hamiltonian should include the kinetic energy operator related with that monomer in the 
Hamiltonian of the cluster and a model potential surface for that monomer. The model 
potential surface can be chosen as desired. However, the efficiency of the resulting basis will 
depend on this choice. 

The primitive basis functions which are used to generate the optimized basis should 
be symmetry adapted by using the operations contained in the pure permutation group 
related with that monomer. According to equation [51 this is certainly sufficient for making 
sequential symmetry adaptation properly. However, in order to guarantee invariance of the 
cluster basis while combining the monomer bases, it is better to follow a more complex 
path for making sequential symmetry adaptation of basis functions because of the reasons 
explained below. 

If the sequential symmetry adaptation procedure is used as it is, then the basis functions 
are going to be symmetry adapted to the irreducible representations of the inversion sub- 
group while forming the symmetry adapted cluster basis functions. At this step, it will be 
necessary to apply inversion operation, E* to the monomer basis functions. When the inver- 
sion operation is applied to a basis function of a monomer, the result will be another basis 
function for the same monomer. However, the resulting function will not be in the same 
basis unless the basis has the inversion symmetry of the cluster. Since the optimized basis 
functions of the monomers are going to be generated with a calculation, the basis functions 
cannot have such a symmetry unless the model potential surface of the cluster has the in- 
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version symmetry of the cluster. A model potential surface for a monomer cannot have such 
a property unless it is imposed to it. Consequently, if the sequential symmetry adaptation 
procedure is used as it is, it will be necessary to deal with a generalized eigenvalue problem 
instead of a standard eigenvalue problem. 

This problem can be overcome by using the properties of direct product groups. The 
inversion subgroup can always be multiplied with another subgroup of the molecular sym- 
metry group with direct product multiplication. If the basis functions that are used to 
generate optimized monomer bases are symmetry adapted to irreducible representations of 
the direct product group obtained from the pure permutation group of the monomer and the 
inversion subgroup of the cluster. Then, they will be symmetry adapted to the irreducible 
representations of both the pure permutation group of the monomer and the inversion sub- 
group of the cluster. In this case, the application of the inversion operation to optimized 
basis functions will not create new functions. In fact, they will be the eigenstates of the 
inversion operation with the eigenvalues ±1. If the same thing is true, for all of the monomer 
bases, then the product basis of the monomer bases will also be eigenstates of the inversion 
operation with the eigenvalues ±1. Consequently, the product basis of the monomer bases 
will be invariant under the effect of the inversion operation so that the symmetry adapta- 
tion of basis functions will not lead to a generalized eigenvalue problem but to a standard 
eigenvalue problem. 

To sum up, in order to generate a small basis for a monomer, the permutation group of 
that monomer, Gk, and the inversion subgroup, e, of the molecular symmetry group of the 
cluster are taken and the direct product group of these two subgroups are formed. Then, 
the eigenstates of the model Hamiltonian is solved for each symmetry separately after the 
basis functions are symmetry adapted. A subset of these eigenstates becomes the contracted 
basis for that monomer. 

B. Generation of Bases for Other Monomers 

The procedure given in section IIII Al can be used for generating a contracted basis for 
each of the monomers. However, this will not be the optimal choice. Because, the molec- 
ular symmetry group of the cluster includes the subgroup which includes the operations 
permuting the identical monomers as a subgroup. Therefore, after the monomer bases are 
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combined together they are going to be symmetry adapted by using the symmetry operations 
contained in this group. The permutation operations in this group will mix the monomer 
bases, such that when these permutation operations are applied to the basis functions of a 
monomer, then the resulting function will be a basis function for another monomer. If this 
resulting basis function is not already available in the contracted basis of that monomer, 
then it will not be orthogonal to the basis functions of that monomer necessarily. Therefore, 
unless there is a relation between the bases of different monomers, symmetry adaptation to 
the group which includes the operations permuting the identical monomers will always be 
a problem. As a result of that, a better way of constructing bases for all of the monomers 
is to find a basis for one of them, and then to generate bases for other monomers from the 
basis of this monomer. 

The obvious choice for generating the bases for other monomers could be just to relabel 
the basis functions of a single monomer for different monomers. However, this will not 
help to get rid of the symmetry adaptation problem posed before. If the monomer basis 
functions are generated with relabeling, the basis functions obtained in this way will be 
symmetry adapted in their own permutation groups since these groups are isomorphic to each 
other. However, if the permutation operations which permute the identical monomers has 
some effect on the coordinates other than relabeling them, the action of these permutation 
operations to basis functions will still result in new basis functions. Although relabeling 
the basis functions of the first monomer to generate bases for the other monomers can be 
a solution for some specific cases, where the action of the permutation operations is just to 
relabel the coordinates (in fact this would be the case if Cartesian coordinates were used), 
it is not applicable in general, and is far from being a general solution. 

The general solution to that symmetry adaptation problem can be found as follows. 
Firstly, let's consider the case of a dimer. If the monomers are labeled as a and b, then the 
group which contains the permutations of identical monomers will be G2 = {E, P a b}, where 
the operation P a b is the permutation operation which permutes the monomers a and b. If 
(p£ is the k th basis functions in the optimized basis of the monomer a and <f>\ is the I th 



function in the optimized basis of the monomer b (which is to be determined), then 



W 



will be one of the product basis functions. In order to make symmetry adaptation of the 
basis functions, it will be necessary to apply the permutation operation P a b to the basis 
functions. The application of the operation P^ will relabel the basis functions, but it will 
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also introduce some changes to them so that the new basis function obtained by application 
of this operation will be 

where <jr k = P a b4>k an< ^ 4>i = Pab<Pi ■ The resulting basis functions are labeled with 's to 
imply that they are not necessarily in the bases of monomers a and b. At this point it can 
be realized that the problem can be get rid of by defining 

<Pf = Pa^- (8) 

Thus, equation © becomes, 

iWi a V, (6> = Pa b {^\Pa b ^)) (9) 

= {Pa^iPabPaJ?) (10) 

= tfV! a) . (11) 

In the equations above P a bPab — E is used that is permuting the two monomers between 
them twice is equivalent to the application of the identity operation, so that it leaves the 
system invariant (In other words, P a b is its own inverse: P ah = P a b)- 

Thus, if the basis functions for the monomer b are generated by using equation (jSJ), then 
the application of the permutation operation P a ^, does not generate new basis functions but 
just carries a basis function in the basis of a monomer to another basis function in the basis 
of the other monomer. 

Although, the discussion above is made for just two monomers, the idea can be extended 
to any bigger cluster. In the case of a trimer, for example, if the monomers are labeled as 
a, b and c, the cyclic group containing the permutations of identical monomers will be the 
group G 3 = {E, P a bc, P ac b}- In this case the basis of the monomer b can be generated by 

$ = Patd? ■ (12) 

The basis of the monomer c can be generated from the basis of monomer b by 

<Pi C) = Pabc^ ■ (13) 

Thus, by repeated application of the generator of the group G3, it is possible to generate 
bases for all of the three monomers from the basis of a single monomer. This method can 
be extended to any bigger cluster. 
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If the group Gi is the cyclic group with the order n m , then a basis for all of the monomers 
can be generated from the basis of a single monomer by repeated application of the generator 
of the group Gi. The reason for the assumption that the group Gi is cyclic should be clear 
at this point. For example, in the case of a trimer if the group containing the permutations 
of the identical monomers would be the group G e = {E, P abc , P ac b, P a b, P ac , Pbc}, there would 
be more than one way of generating bases for other monomers. For example, both of the 
operations P a b c and P a b can be used to generate a basis for the monomer b from the basis of 
the monomer a. If the operation P a b c (P a b) is used; then, the application of the operation P a b 
(P a bc) ma y still create new basis functions. In such a situation, it is impossible to guarantee 
the invariance of the basis. 

Before closing this section, it should also be noted that a basis, which is generated by using 
the generator of the group which includes the operations that permute identical monomers, 
will have the same orthogonality relations with the original basis. For example, if one has an 
orthonormal basis for the monomer a, then the basis of the monomer b generated by using 
equation (JE]) will have 

(<pf\<P?) = (PaA^Pab^) 
= i^lPlbPab^) 

= (4°W 0) > 

= hi (14) 

Besides, the basis functions of the monomer b will be eigenstates of the model Hamiltonian 
of the monomer b which is generated in the way that the eigenstates of the monomer b is 
generated. Thus, if H® is the model Hamiltonian of the monomer a, and cp£ is the k th 
eigenstate of this model Hamiltonian with the eigenvalue ej^ such that 

H>t ] = 4 Q) 0i a) ; (15) 

then, by applying the permutation operation to both sides of the equation above, 

P ab Hyt ] = Pabei a) <f>t\ (16) 

and inserting the identity operation P ab P a b = E between Hamiltonian and the basis function 
one gets 

PabHAPab^ = e { ^ P ab ^ . (17) 
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Thus, 

H^f = e[ a Ut\ (18) 

where equation (jSJ) is used , and H® is defined as 

H° b = P ab H° a Pl b . (19) 

Therefore, from equation ffl8|) . it follows that (fi k is an eigenstate of the model Hamiltonian 
H® with the eigenvalue e k = e k , such that 

^Vi 6) =4 6) 0i 6) - (20) 

C. Combining Monomer Bases 

After a basis for one of the monomers is generated by solving for the eigenstates of a 
model Hamiltonian, and the bases for the other monomers are generated from the basis of 
this monomer by using the generator of the group containing the permutations of identical 
monomers; a basis for calculating the spectra of the cluster can be generated by forming 
the tensor product of the monomer bases and combining them with a primitive basis for the 
inter-monomer coordinates. Thus, if (f>nj is the n th basis function for the monomer k and xi 
is the I th basis function related with inter-monomer coordinates, the eigenstates of the full 
problem can be expanded in this product basis as 

n m 

%= E ^r ,n2 '-' nnm, ^nC' ( 2l ) 

ni,ri2,...,n nm ,l k=l 

where (7™ 1 > n2 >-> n ">«> j g ^q expansion coefficient. 

As discussed in the previous sections, the product basis of the monomer bases will be 
invariant under the effect of any permutation inversion operation. Since a primitive basis 
is used for the inter-monomer coordinates, it can be chosen to be invariant under the effect 
of any permutation inversion operation. Consequently, the eigenstates of the cluster can 
be obtained as linear combinations of the basis functions of this basis even if the basis is 
not symmetry adapted for calculations. However, symmetry adaptation is always useful for 
reducing the size of the basis. 

The obvious way to make the symmetry adaptation of the basis functions to an irreducible 
representation of the molecular symmetry group of a cluster is to apply the projection 
operator of that irreducible representation to the basis functions. 
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Since the group Gk®£ which is used in monomer calculations contains the group Gk as an 
invariant subgroup, the monomer basis functions will already be symmetry adapted in their 
own pure permutation groups G k . Therefore, symmetry adaptation of the basis functions 
can be done in two steps after finding the correlations between the irreducible representations 
of the molecular symmetry group and the irreducible representations of the pure permutation 
groups. For example, considering the bases of the monomer a if the irreducible representation 
r a of the molecular symmetry group correlates to the irreducible representation Tj of the 
group G k , then only the bases which are symmetry adapted to the irreducible representation 
Tj should be used for the calculations. The bases which are symmetry adapted to the other 
irreducible representations will be annihilated by the application of the projection operator 

of the group T a , since this projection operator contains the projection operator P J . of the 

G k 

group G k . On the other hand, the application of the projection operator P \ a) will leave 
the basis functions of the monomer a which are symmetry adapted to the Tj irreducible 
representation of the group G k invariant. Thus, if the correlations between the irreducible 
representations of the molecular symmetry group of the cluster and the pure permutation 
groups of the monomers are found and the appropriate bases are chosen, then the symmetry 
adaptation of the basis functions can be done in two steps by symmetry adapting the basis 
functions further to the groups Gi which is the group containing the permutations of the 
identical monomers and the inversion group e. 

Since monomer basis functions are symmetry adapted to the group Gk <8> e, there will 
be two bases which are symmetry adapted to two different irreducible representations of 
the group Gk <S> s and at the same time which are symmetry adapted to the irreducible 
representation Tj of the group Gk- One of these bases will be symmetry adapted to the 
irreducible representation Yi ® G = Ti g , which will have even parity, and one of them will 
be symmetry adapted to the irreducible representation F{ (g) U — Ti U (see table [III for the 
irreducible representations G and U of the group e). Thus, if it is necessary to have a basis 
for the monomer a which is symmetry adapted to the irreducible representation I\, then the 
basis which should be used for the monomer a will be 

*- ig tt7 *■ iui 

where the labels of the irreducible representations are used to mean any basis function 
belonging to that symmetry. The bases of all of the monomers can be found similarly. Thus, 
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when the correlations are found and the product basis is formed there will be 2™ m different 
product bases which differ from each other by the symmetries of monomer functions. Since 
inter-monomer coordinates are usually invariant under the effect of the inversion operation, 
half of these terms will have even parity and half of them will have odd parity. Thus, when 
the symmetry adaptation of the basis functions to the inversion group is done half of these 
terms will be annihilated and half of them will leave invariant. After finding the correlations, 
a fully symmetry adapted basis can be generated by symmetry adapting them by using the 
symmetry operations of the pure permutation group Gi. Thus, the symmetry adapted basis 
functions can be obtained as 

1 Z| g&Gt fc=l 

where T a is an irreducible representation of the molecular symmetry group, and x rQ! [<7] 
is the character of the operation g represented with the operation O g in the irreducible 
representation T a . 

D. Solution of The Full Problem 

In order to find the eigenvalues of the Hamiltonian, the matrix elements of the matrix 
representing the Hamiltonian should be calculated. They can be evaluated easily if the 
Hamiltonian is partitioned as 

Tim 

H = J2 H°k + AT + A ^> (23) 

fc=i 

where H^s are the model Hamiltonians for the monomers, AT is the kinetic energy terms 

which are not included in the model Hamiltonians and AV = V — Y^t^i ^4° ^ s ^ ne difference 

between the potential surface of the full problem and the sum of the model potential surfaces 

used in the model Hamiltonians to determine the bases for the monomers. 

The basis functions of the full problem become eigenstates of the zeroth order Hamiltonian 

for the full problem, such that by defining 

H° = J2 H ° k , (24) 

fc=i 

and expanding the wave function in the product basis 

' *■ m 

^n u n 2 ...,n nm ,l = XlH^i (25) 



fe=l 



15 



the following eigenvalue relation is obtained 

/ j e n k I Yn 1 ,n 2 ,—,n nm ,l- (,2b) 

fe=l / 

Therefore, the matrix elements of the Hamiltonian in the basis of the full problem will be 
given by 



i^n' 1 ,n' 2 ,...,n' nm ,v\ H \ % l ) n 1 ,n 2 ,...,n nm ,l) = / ] e n k .hi' [j S n r n' r 

k=l r=l 

+ (W 1 ^,...,< m ,H A " T + £V\if> ni)nar .. )nnm>l ). (27) 



Thus, in order to calculate the matrix elements of the Hamiltonian, it is necessary to evaluate 
matrix elements of the AT and AV terms, in the basis of the full problem. These terms can 
be evaluated in the primitive bases of monomers and in the primitive basis of inter-monomer 
coordinates, then they can be transformed to the contracted basis of the cluster. If the terms 
AT and AV are small; then these terms can be considered as a small perturbation and the 
basis functions will resemble to the eigenstates of the actually problem. In such a case, 
convergence of the results can be obtained by using a small number of contracted basis 
functions. However, this may not be the case for many problems. 

IV. FURTHER DETAILS ABOUT THE MBR METHOD 
A. The Case of Non-identical Monomers 

While discussing the method in section IHIl it was assumed that the cluster consist of 
identical monomers. However, this assumption can be relaxed easily. 

If there are non-identical monomers in the cluster, then a monomer basis can be generated 
for each type of monomer separately in the way it is discussed in section IIII A[ In such a 
situation, for each type of monomer, there will be subgroups of the molecular symmetry 
group that contain permutation operations which permute identical monomers. The group 
Gi which is the group that contains the permutations of identical monomers will be the 
direct product of these subgroups. 

After a basis for a single monomer of one type of monomers is generated, the procedure 
of section IIII Bl can be used to generate bases for other monomers which are identical with 
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that monomer. The same thing can be done for each type of monomers, so that a basis for 
each monomer in the cluster is generated. 

Once a basis for each monomer is generated, the procedures of sections IIIICI and IIIIDI 
can be used to combine the bases and to solve the problem. These steps has nothing to do 
with whether the monomers are identical or not. 

B. Use of The Method with Pseudospectral Methods 

Pseudospectral methods are used frequently in quantum mechanical calculations. In this 
type of methods, two different bases, which are usually isomorphic to each other, are used 
for the same degrees of freedom. Then, the different parts of the Hamiltonian are evaluated 
in the basis whichever is convenient for that part of the Hamiltonian. 

If pseudospectral methods are used for evaluating the matrix elements of the Hamiltonian 
operator, then it becomes necessary to make transformations from one basis to another. This 
is achieved by using transformation matrices. For the monomer a, if <pf with i = 1, 2, . . . , N 
are basis functions in one of the bases and Of 1 with i = 1,2, ... ,N are the basis functions 
in the second basis, then the relation between the basis functions of the two different bases 
will be 

j=N 

4 a) = E T ^ a) ( 28 ) 

where T^ are the matrix elements of the transformation matrix T. 

In section IIHB[ it was seen that the proper way of generating bases for other monomers 
is to use the generator of the group which contains the operations that permute identical 
monomers. If the permutation operation P a j, is applied to the both sides of the equation 
above, by using equation flB]) and the fact that matrix elements are just constants, the 
equation becomes 

j=N 

Since the basis functions (j)\ and <pf are related to each other by equation OH]) consistency 
requires that the same thing should be true for the basis functions 6\ and 9\ , that is 

a( b ) _ r> n(a) 



PabOr- (30) 
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Thus, equation (1291) becomes 

j=N 

4 fe) = E T ^ 6) - ( 31 ) 

3=1 

Equations (1281) and (131]) show that if the basis functions of the monomers are related to each 
other with equation ([H]); then, the matrix elements of the transformation matrix for the bases 
of the monomer b will be the same with the matrix elements of the transformation matrix 
for the bases of the monomer a, so that both of the transformations can be done with the 
same transformation matrix. Consequently, while doing computations, it is not necessary to 
store separate transformation matrices for each monomer. The same transformation matrix 
can be used for all of the monomers of the same type. Since the transformation matrices 
usually occupy large memories especially if multidimensional coupled bases are used, this 
fact is very useful for reducing the memory cost of computations. 

V. QUALITATIVE MODELS 

In the following two subsections, two qualitative models will be given for the possible 
applications of the MBR method. The main discussion will be about the symmetries of 
wave functions. An implementation of the qualitative model given for water dimer will be 
made in a following paper for calculating its vibration-rotation-tunneling spectra. 

A. Water Dimer 

A group theoretical treatment of water dimer is done by Dyke and co-workers for ex- 
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plaining the microwave data |9| . The molecular symmetry group of this dimer is the group 
Gi6 which is isomorphic to the D^ point group. If the oxygen atoms in the molecule are 
labeled as a and b; the hydrogen atoms bonded to the oxygen a are labeled as 1 and 2; and 
the hydrogen atoms bonded to the oxygen b are labeled as 3 and 4; then, this group can be 
written as 11] 

G 16 = ((Gt } ® Gf ) ©G ( 2 ab) ) (8) e, (32) 

where the monomer permutation groups are 

G$° = {£,(12)}, (33) 

Gf = {£,(34)}, (34) 
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TABLE I: Character table of the permutation group G 2 . 



G 2 E (12) 


A 1 1 
B 1 -1 



TABLE II: Character table of the inversion group e. 



e E E* 
G 1 1 
U 1 -1 



the group containing the operations that permute the identical monomers is 

Gf ) = {E,(a6)(13)(24)} ! (35) 

and the inversion group is 

e = {E,E*}. (36) 

Character table of the group G 2 is given in table [I] Character table of the group G 2 
can be obtained from the character table of the group G 2 by replacing the operation (12) 
with the operation (34), since the two groups are isomorphic to each other. Similarly, the 
character table of the group G 2 can be obtained from the character table of the group 
G 2 by replacing the operation (12) with the operation (a6)(13)(24). Character table of the 
group e is given in table [TT1 

The monomer calculations should be done by symmetry adapting the basis functions to 
the irreducible representations of the group which is obtained as a direct product of the pure 
permutation group of the monomer and the inversion subgroup of the cluster. In the case 
of water dimer, if the monomer is chosen as the monomer a, this means that the monomer 
calculations should be done by symmetry adapting the basis functions to the irreducible 
representations of the group 

C 2v (M) = G { 2 a) ®e = {E, (12), E*, (12)*}. (37) 

Character table of the group C 2v (M) is given in table IIHI 
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TABLE III: Character table of the C 2v (M) permutation inversion group. This group is the direct 
product of the groups G 2 whose character table is given in table U and the group e whose character 
table is given in table HH In the table, T = x ® y means that the irreducible representation T of 
the group C 2v (M) is obtained as direct product of the irreducible representation x of the group 
G 2 and the irreducible representation y of the group e. 



C 2v {M) = G 2 ®e 


E 


(12) E* (12)* 


A x = A ® G 


1 


1 1 1 


A 2 = A®U 


1 


1 -1 -1 


B 1 = B®U 


1 


-1 -1 1 


B 2 = B®G 


1 


-1 1 -1 



After the monomer calculations are done with the basis which is symmetry adapted to the 
irreducible representations of the group C 2v (M), the basis of the monomer b can be generated 
by using the generator of the group G 2 which is (afe)(13)(24). Then, the question becomes 
how to combine these bases for the solution of the full problem. Character table of the 
group Gi6, which is the molecular symmetry group of water dimer, is given in table |V] The 
correlations between the irreducible representations of the group G\§ and the irreducible 
representations of its subgroups are given in table IVll 

The monomer bases will already be symmetry adapted to the irreducible representations 
of the groups G 2 and G 2 ■ Therefore, the product basis of monomer bases will be symmetry 
adapted in the group G4 = G 2 <8> G 2 , whose character table is given in table IIVI 

Considering the A^ representation of the group Giq, the correlation table given in table lVIl 
says that this level correlates to I\ irreducible representation of the group G4, A irreducible 
representation of the group G 2 and G irreducible representation of the group e. Therefore, 
if there is a basis for the monomers such that its basis functions are already symmetry 
adapted to the I\ representation of the group G4, then these basis functions can be symmetry 
adapted to the irreducible representations of the group G\q by applying the projection 
operators of the A representation of the group G 2 and the G representation of the group 
e. 

The I\ irreducible representation of the group G4 is the direct product of the A repre- 
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TABLE IV: Character table of the group G4. In the table, r = x ® y means that the irreducible 
representation F of the group G4 is obtained as a direct product of the irreducible representation 
x of the group G 2 an d the irreducible representation y of the group G 2 ■ 



G 4 = G { 2 a) ® Gf 


E 


(12) (34) (12)(34) 


Ti = A® A 


1 


1 1 1 


T 2 = A®B 


1 


1 -1 -1 


f 3 = b®a 


1 


-1 1 -1 


F A = B®B 


1 


-1 -1 1 



TABLE V: Character table of the G\% PI group which is the molecular symmetry group of water 
dimer. This group is isomorphic to the D A h point group. This character table is taken from a 
paper by Dyke |9(. 





(12) (a6)(13)(24) (o6)(1324) 






(12)* 


(a6)(13)(24)* 


(a6)(1324)* 




Gi6 E 


(34) (a6)(14)(23) (o6)(1423) (12) (34) 


E* 


(34)* 


(a6)(14)(23)* 


(a6)(1423)* 


(12)(34)* 


At 1 


1 


1 1 


1 


1 


1 


1 


1 


1 


A+ 1 


-1 


-1 1 


1 


1 


-1 


-1 


1 


1 


Bt 1 


1 


-1 -1 


1 


1 


1 


-1 


-1 


1 


B+ 1 


-1 


1 -1 


1 


1 


-1 


1 


-1 


1 


E + 2 








-2 


2 











-2 


A', 1 


1 


1 1 


1 


-1 


-1 


-1 


-1 


-1 


A 2 1 


-1 


-1 1 


1 


-1 


1 


1 


-1 


-1 


sr 1 


1 


-1 -1 


1 


-1 


-1 


1 


1 


-1 


B 2 ~ 1 


-1 


1 -1 


1 


-1 


1 


-1 


1 


-1 


E~ 2 








-2 


-2 
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sentation of the group G 2 and the A representation of the group G 2 ■ Therefore, the direct 
product of the monomer bases which are symmetry adapted to the A representations in their 
own permutation groups will be symmetry adapted to the Ti irreducible representation of 
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TABLE VI: Correlations between the irreducible representations of the group G\q and its sub- 
groups. 



16 ^4 ^2 ^ 



Af T ± 


A G 


At r 4 


5 G 


b+ ri 


B G 


b+ r 4 


A G 


e+ r 2 e r 3 


A®B G 


^r Tx 


A U 


^ r 4 


B U 


b x - i\ 


B U 


B 2 r 4 


A U 



E~ To © T 3 A®B U 



the group G4. 

The monomer calculations are done with the group C2„(M). In this group both of 
the irreducible representations A 1 and A 2 contains the A representation of the group G 2 > 
therefore both of the bases which are symmetry adapted to the irreducible representations 
A\ and A2 of the group C2 V {M) can be used for the monomer a, so that the basis of the 
monomer a should be Ai © A 2 , where the labels of the irreducible representations are used 
to imply basis functions belonging to that symmetry. The same thing will be true also for 
the monomer b. Thus, there will be four different product bases which will be symmetry 
adapted to the irreducible representation 1^ of the group G4 which are: A\ <g> A±, A 1 <g> A 2 , 
A 2 ®A U A 2 ® A 2 . 

The character of the operation E* is 1 in the irreducible representation A\ and —1 in the 
irreducible representation A 2 . Therefore, the basis functions included in the product bases 
A\ ®A\ and A 2 <g> A 2 will not change sign upon the operation of the inversion operation E*, 
so that they will have even parity. On the other hand the basis functions included in the 
product bases A\ ® A 2 and A 2 ® A\ will change sign upon the operation of the inversion 
operation E*, so that they will have odd parity. In order to fully symmetry adapt the these 
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TABLE VII: This table shows which monomer bases should be combined for obtaining bases for 
water dimer calculations with the group Gxe- In the table, labels of the irreducible representations 
are used to imply basis functions belonging to that symmetry. For an explanation of how to obtain 
mutually orthogonal basis for the doubly degenerate levels (i.e. Ej, Ey) see reference 
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Gxe Bases Gxe Bases 

Af (Ax © A t ) © (A 2 x A 2 ) A{ (A x © A 2 ) © (A 2 © A x ) 
A+ (B x © Si) © (B 2 © B 2 ) A^ (B x © B 2 ) © (B 2 © Bi) 
S+ (yli Ai) © (A 2 © A 2 ) Sf (A t © A 2 ) © (A 2 © 40 
S+ (Si © B x ) © (S 2 © S 2 ) J3J (Si © S 2 ) © (S 2 ® Si) 
S+ (Ai © S 2 ) © (4 © Bx) E~ (Ax © Si) © (A 2 © S 2 ) 
K+ (B 2 © Ai) © (Si © A 2 ) E- (Bx © Ax) © (S 2 © A 2 ) 

basis functions to the Af irreducible representation of the group Gxe, it is necessary to 
apply the projection operator of the G irreducible representation of the group e, which is 
P e = (E + E*)/2. Since the basis functions in the product bases Ax <8> Ax and A 2 ® A 2 does 
not change sign upon the operation of the inversion operator, these basis functions will be 
invariant under the effect of the projection operator P e G . On the other hand, since the basis 
functions in the product bases Ax £5> A 2 and A 2 <g> Ax change sign upon the operation of the 
inversion operation, they will be annihilated when the projection operator P G is applied. 
Therefore, among the four product bases which are symmetry adapted to the I\ irreducible 
representation of the group r 1; only the bases Ax ® Ax and A 2 ® A 2 can be used for the 
calculations of the Af representation of the group Gxe- 

Thus, the way to generate a basis for the calculations of the Af representations of the 
group Gxe is to take the product bases Ax <8> Ax and A 2 ® A 2 and symmetry adapt them 
to the irreducible representation Af of the group Gxe by applying the projection operator 



P\ ab ) of the group G 2 . Application of the projection operator P G will already leave the 



i(ab) 
r 2 • 

basis functions invariant. 

Similarly, the way to combine the monomer bases for the other irreducible representations 
of the group Gxe can be found. The results are summarized in table IVIII 
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TABLE VIII: Character table of the group G24. This table is taken from a paper by van der 
Avoird et. al. [15J]. In the table, representations given as A2 g , A% g , A2 U and ^3^ in that article 
are combined to the doubly degenerate representations E g = A2 g © A% g and E u = A2 U © A% u . In 
contrast to the reference 15J], operations in the classes are shown explicitly. 



{acb) (164253) (abc) (135264) (acb) (153) (264) (abc) (135) (246) 

(ac6)(153264) (a6c)(146235) (acb) (164) (253) (a&c) (146) (235) (12) (12)(34) 

(ac6)(154263) (a6c)(145236) (acb) (154) (263) (abc) (145) (236) (34) (12)(56) 

E (ac6)(163254) (a6c)(136245) (acb) (163) (254) (abc) (136) (245) (56) (34)(56) (12)(34)(56) 



Mg 1 


1 


1 


1 


1 


1 


1 


1 


Eg 2 


-1 


-1 


-1 


-1 


2 


2 


2 


T g 3 














-1 


-1 


3 


Am 1 


-1 


-1 


1 


1 


-1 


1 


-1 


E u 2 


1 


1 


-1 


-1 


-2 


2 


-2 


T u 3 














1 


-1 


-3 



B. Water Trimer 



A group theoretical treatment of water trimer is done by van der Avoird et. al. 15[. The 
molecular symmetry group of water trimer is the group G^. This group is the direct product 
of the inversion subgroup and the pure permutation subgroup G24, such that G48 = G24 <E>£- 
Character table of the group G24 is given in table IVIIII In the table a, b, c are the labels 
for the oxygen atoms. The hydrogen atoms bonded to oxygen a are labeled as 1 and 2, the 
hydrogen atoms bonded to oxygen b are labeled as 3 and 4, and the hydrogen atoms bonded 
to oxygen c are labeled as 5 and 6. 

If follows from equation fl3]), the molecular symmetry group of water trimer can be written 
in terms of its subgroups as 11] 



G 



18 



g { : ] 



Gf 



y(c) 



GV ©G 3 ®e 



(38) 
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where the pure permutation groups of the monomers are 

G 2 a) = {£1,(12)}, (39) 

Gf = {£,(34)}, (40) 

Gf = {£,(56)}, (41) 

and the group containing the operations permuting the monomers is 

G 3 = {£, (a&c) (135) (246), (oc6)(153)(264)}. (42) 

Character table of the group G% is given in table IIXI 

TABLE IX: Character table of the group G3. In the table, w = exp(2i7r/3). 



G 3 


E 


(06c) (135) (246) (acb) (153) (264) 


A 


1 


1 


1 


E, r 


1 


w 


w* 


Ey 


1 


w* 


w 



As long as the symmetry adaptation of the monomer calculations are concerned, there 
is no difference between water dimer and water trimer. The monomer calculations should 
be done by symmetry adapting the basis functions to the irreducible representations of the 
group C2v(M) = G 2 ® e. After the monomer calculations are done, bases of monomer 
b can be generated from the bases of monomer a by using the permutation operation 
(abc) (135) (246). The bases of the monomer c can be generated either from the bases of 
the monomer a by using the permutation operation (ad)) (153) (264) or from the bases of the 
monomer b by using the permutation operation (abc) (135) (246). 

Since the monomer bases functions are already symmetry adapted in their own permu- 
tation groups, the product basis of the monomer bases will be symmetry adapted in the 
group G% = G 2 ® G 2 ® G 2 . Character table of the group Gg is given in table |2D The 
correlations between the irreducible representations of the group G 48 and its subgroups G$, 
G3 and e are given in table IXIl 

In order to illustrate how the monomer bases can be combined for the solution of the full 
problem consider the £ p h representation. This representation correlates to the 1^ representa- 
tion of the group Gg. Since the irreducible representation r\ is obtained as a direct product 
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TABLE X: Character table of the group Gg. Correlations with the irreducible representations 
of the subgroups are indicated such that in the table r = x <8> y ® z means that the irreducible 
representation T of the group Gs is obtained as a direct product of the irreducible representations 
x, y and z of the groups G 2 , G 2 and G 2 , respectively. 



G 8 = G 2 a) ®G i 2 b) ®G 2 c) 


E 


(12) (34) (56) (12)(34) (12)(56) (34) (56) (12)(34)(56) 


Tx = A <8 A (g) A 


1 


1 


1 


1 


1 


1 


1 1 


T 2 = A <8 A (g) B 


1 


1 


1 


-1 


1 


-1 


-1 -1 


T 3 = A ig) S (8) A 


1 


1 


-1 


1 


-1 


1 


-1 -1 


T 4 = B(g>A®A 


1 


-1 


1 


1 


-1 


-1 


1 -1 


T 5 = A®B®B 


1 


1 


-1 


-1 


-1 


-1 


1 1 


T 6 = B®A®B 


1 


-1 


1 


-1 


-1 


1 


-1 1 


T 7 = B®B®A 


1 


-1 


-1 


1 


1 


-1 


-1 1 


r 8 = B®B®5 


1 


-1 


-1 


-1 


1 


1 


1 -1 



of the A representations of the pure permutation groups of the monomers (see table |X"I) . the 
product basis of the monomer bases which are symmetry adapted to the A representation 
should be used for forming the monomer bases for the solution of the full problem. Since 
the monomer calculations are done with the group C 2v (M), both of the A\ and A 2 repre- 
sentations are symmetry adapted to the A representation of the pure permutation group of 
the monomers. Therefore, the monomer bases that should be used in the calculations is the 
sum of A\ and A 2 bases which is Ai © A 2 , and the product basis of the monomer bases is 
going to be 

{A 1 © A 2 ) © {At © A 2 ) © (Ai © A 2 ). 

When the multiplications are done, there will be eight terms. Half of them will have even 
parity and half of them will have odd parity. The functions having odd parity will be 
annihilated upon the operation of the projection operator P e G of the inversion group. On 
the other hand the functions having even parity will be invariant under the effect of the 
projection operator P e G . Thus, only the terms which have the right parity will lead to basis 
for the full problem. These bases are {A x © A x © A x ) © (A x © A 2 © A 2 ) © (A 2 © A x © A 2 ) © 
(A 2 © A 2 © A±). After these bases are combined with a basis related with inter-monomer 
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TABLE XI: Correlations between the irreducible representations of the group G 48 and its sub- 
groups. 



G48 G 8 


G 3 


e 


At 9 Ti 


A 


G 


e+ 2 * ri 


E x @E y 


G 


t+ r 5 e r 6 e r 7 


A@E X ® Ey 


G 


Af n r 8 


A 


G 


e+ 2 * r 8 


E X ffi Ey 


G 


t+ r 2 e r 3 e r 4 


A®E x ®E y 


G 


K 9 ?i 


A 


U 


E- 2 * Tx 


E x @E y 


u 


t- r 5 er 6 er 7 


A@E x ®E y 


u 


^i~« r 8 


A 


u 


e- 2 * r 8 


E X © Ey 


u 


t- r 2 er 3 er 4 


A@E x ®E y 


u 



coordinates, they can be symmetry adapted to the E+ irreducible representation of the 
group G 48 by application of the projection operator P G * + P G ^- 

Similarly, the monomer bases which should be used to form a product basis for each 
irreducible representation of the group G 48 can be found. The results are given in table Kill 

VI. CONCLUSIONS 

By using the sequential symmetry adaptation procedure, and the physically meaningful 
partitioning of the molecular symmetry group of molecular clusters, a new method for calcu- 
lating the VRT spectra of molecular clusters named Monomer Basis Representation (MBR) 
method was developed in section IIIII In the MBR method, calculations starts with a single 
monomer with the purpose of obtaining an optimized basis for that monomer as a linear 
combination of some primitive basis functions. Then, an optimized basis for each identical 
monomer is generated from the optimized basis of this monomer. By using the optimized 
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TABLE XII: This table shows which monomer bases should be combined for obtaining bases for 
the calculations of the water trimer with the group G48. In the table labels of the irreducible 
representations are used to imply the monomer basis functions belonging to that symmetry. For 
an explanation of how to obtain mutually orthogonal bases for the triply degenerate levels (i.e. 
T 9 t; T gyi T gz), see reference [u|. 



G48 Bases 



A + 


[Ai ® Ai® Ai 


) © (Ax ® A 2 ® A 2 


) © (A 2 ® Ax ® A 2 


)ffi(A 2 ® A 2 ® Ai) 


^ 


A\ ®A x ®Ax 


© (Ai (8 A 2 (8) A 2j 


© (A 2 ® Ai ® a 2 ; 


©(A 2 ® A 2 ® Ai) 


T+ ( 


Ai®B 2 ® B 2 ] 


© (Ax ® Ei ® Bx] 


©(A 2 ®E 2 ®Ei; 


© (A 2 ® Ei ® E 2 ) 


T+ ( 
gy y 


B 2 ®B 2 ®Ax] 


© (Bx ® Bx ® Ai; 


© (E 2 ® Ei ® A 2 ] 


© (Ei ® E 2 ® A 2 ) 


T+ ( 
gz 1 


B 2 ®A 1 ® B 2 ] 


(Si ® Ai ® E^ 


© (Ei ® A 2 ® e 2 ; 


ffi(E 2 ® A 2 ®Ei) 


A + 


[B 2 ®B 2 ®B 2 


) e (s 2 ® -Bi ® e x 


) © (Ei ® E 2 ® Ei 


) © (Ei ® Ei ® E 2 ) 


E+( 


B 2 ®B 2 ® B 2 ] 


© (E 2 ® -Si <8 Bx] 


© (Ei ® e 2 ® e x ; 


© (Ei ® Ei ® E 2 ) 


T + 
ux 


[Ai ®A 1 ®B 2 


) © (Ai ® a 2 ® 5i 


) © (A 2 ® Ax ® E 2 


© (A 2 ® A 2 ® E 2 ) 


uy 


A 1 ®B 2 (g) A x 


e (A 2 ® £1 ® Ai 


© (Ai ® E 2 ® a 2 ; 


© (A 2 ® E 2 ® A 2 ) 




B 2 ® A x ® Ax 


© (Ei ® Ai ® A 2 


©(e 2 ® a 2 ® a x ; 


© (E 2 ® A 2 ® A 2 ) 


K 


[A 2 ® A 2 ® A 2 


) © (Ai ® Ax ® A 2 


) © (Ai ® A 2 ® Ai 


)ffi(A 2 ® Ai® Ai) 


E g 


A 2 ® A 2 ® A 2 


© (Ax ® a x ® a 2 ; 


ffi(Ai® A 2 ® a x ; 


©(A 2 ® Ai® Ai) 


T~ ( 
gx [ 


A 2 ®Bx® B ±i 


© (Ai ® e 2 ® e x ; 


© (Ai ® Ei ® B 2 ] 


© (A 2 ® E 2 ® E 2 ) 


T~ ( 

gy y 


Bi®Bi® A 2/ 


© (E 2 ® -Bi ® Ax] 


ffi(Ei®E 2 ® a x ; 


© (E 2 ® E 2 ® A 2 ) 


T~ ( 

gz \ 


Bx®A 2 ® Bx) 


© (Si ® Ax ® e 2 ; 


ffi(E 2 ® Ai®Ei; 


© (E 2 ® A 2 ® E 2 ) 


A~ lu 


[Bx ®Bx®Bx 


) © (E 2 ® E 2 ® Bx 


) © (E 2 ® Ei ® E 2 


) © (Ei ® E 2 ® E 2 ) 


K( 


Ei ® Ei ® Bx] 


®(B 2 ®B 2 ®Bx] 


© (E 2 ® Ei ® B 2 ] 
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bases of the monomers, a basis is generated for the solution of the full problem, and the 
VRT spectra of the cluster is obtained by using this basis. Since an optimized basis is used 
for each monomer which has a much smaller size than the primitive basis from which the 
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optimized bases are generated, the MBR method leads to an exponential decrease in the 
size of the basis that is required for the convergence of the results. 

The generation of an optimized basis for each monomer as a linear combination of some 
primitive basis functions led to two problems related with the symmetry adaptation. When 
a basis for the solution of the full problem is obtained as a direct product of the optimized 
monomer bases. These basis functions should still be symmetry adapted to the inversion 
group and also to the group which contains the permutations of the monomers. In order 
to symmetry adapt basis functions to these groups, it is necessary to apply the symmetry 
operations contained in these groups to the basis functions. When these operations are 
applied to the basis functions, the resulting functions should already be in the basis in 
order to have an orthogonal basis. However, when the monomer basis functions are linear 
combinations of some primitive basis functions, the functions obtained as a result of the 
operations of the permutation inversion operations will not be necessarily in the same basis. 

It has been suggested that the invariance of the basis functions under the effect of the 
inversion operation can be achieved by symmetry adapting primitive monomer basis func- 
tions to the inversion subgroup while generating the optimized basis for the monomers. If 
the monomer basis functions are already symmetry adapted to the inversion subgroup, then 
they will be eigenstates of the inversion operator with the eigenvalues ±1. Consequently, 
the basis functions which are obtained as a direct product of the monomer basis functions 
will be the eigenstates of the inversion operator with the eigenvalues ±1. Thus, the basis 
which is obtained as a direct product of the optimized monomer bases becomes invariant 
under the effect of the inversion operation. 

It has also been shown that the product basis of the optimized monomer bases can be 
made invariant under the effect of the operations permuting identical monomers by finding 
an optimized basis for a single monomer and then generating bases for other monomers from 
the basis of that monomer. If an optimized basis for each monomer is generated separately 
independent of other monomer bases, then there is no way to guarantee the invariance of the 
basis. On the other hand, if all of the monomer bases are generated from a single monomer 
basis, then it is possible to correlate the bases to each other so that the product basis of 
optimized monomer bases becomes invariant. It has been shown that a basis for all of the 
monomers can be generated from the basis on a single monomer by repeated application of 
the generator of the group Gi which is the group containing the operations that permute 
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identical monomers. This way of generating bases for all of the monomers made it necessary 
to assume that the order of the group G\ is equal to the number of the monomers: G\ = n m . 

While developing the method, the primitive bases of the monomers were never referenced. 
Thus, the nature of the primitive monomer bases does not matter. They can be product bases 
of one dimensional bases, or they can be multidimensional coupled bases. The symmetry 
adaptation procedure given here will work regardless of the nature of the primitive bases 
used in the monomer calculations. 

Another point which is not discussed in detail in the development of the method was 
what should be the model potential surface while generating a monomer basis. This is 
just because it is impossible to suggest a perfect model potential surface regardless of the 
particular problem being studied. The model potential surface should be chosen according 
to the problem at hand. Nevertheless it should be noted that the more the model potential 
resembles to the actual problem, the more efficient the monomer bases will be. 

Decomposition Of Projection Operators 

In order to find a relation between the projection operators of a product group and the 
projection operators of its subgroups, it is necessary to find a relation between the characters 
of the elements of the product group and the characters of the elements of the subgroups. 
Such a relation can be derived from the following character equation which holds in any 
irreducible representation of any group [161 ]: 

NixWNtxiCj) = dJ2c ijk N kX (C k ). (43) 

k 

In the equation above, Ci, Cj and C k are classes of the group; Ni, Nj and N k are the number 
of elements in these classes; and the coefficients cy* are defined by the class multiplication 
equation CiCj = ^2 lk Cij k C k \ and d is the dimension of the irreducible representation. 

Consider a group G that can be written as a semi-direct product of two of its subgroups 
H and K such that G = H(§)K. If h is an element of group H that is in class Gi of the 
group G and k is an element of the group K that is in class Cj of the group G, and g = hk 
is an element of the group G that is in class C m ; then provided that the class multiplication 
constants satisfy the equation 

Cijk = rS km , (44) 
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where r is an integer; equation ( 143|) reduces to 

X (h)x(k) = d X (g), (45) 

and equation ([5]) follows. Consequently, equation (144)) is the sufficient condition that the 
sequential symmetry adaptation will work for any irreducible representation of the product 
group. This equation seems to be satisfied in many cases for the physically meaningful 
partitioning of the molecular molecular symmetry group of molecular clusters given in equa- 
tion ([3]). The examples include molecular symmetry groups of the clusters (H 2 0) 2 , (C0 2 )2, 
(H20)3, (i?20)2-D2O. However, although it has been argued before [ll| that the equation 
(jUJ) holds for any irreducible representation of any semi-direct product group, this is not 
the case. For example, if pairwise permutations of monomers would be a feasible symmetry 
operation for water trimer so that the molecular symmetry group would be Gqq instead of 
G48, then equation (|4~4"I) would not hold for the semi-direct product multiplication defined 
by equation ([3]). 

It should also be noted that equation (H5|) holds for any one dimensional representa- 
tion. This follows from the fact that one dimensional representations are representations by 
nonzero complex numbers and a representation should satisfy equation fj4"5"]) by definition 
of representation since the characters of the elements are just the complex numbers rep- 
resenting them. Consequently, even if the condition given in equation ( 1441) does not hold, 
sequential symmetry adaptation can still be used for one dimensional representations. 
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